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^ ■ Abstract 

C^^ . Collective excitations spectrum of Dirac electrons in mono and bilayer graphene in the presence 

of a uniform magnetic field is investigated. Analytical results for inter-Landau band plasmon 
j^ ' spectrum within the self-consistent-field approach are obtained. SdH-type oscillations that are a 

' ■ monotonic function of the magnetic field are observed in the plasmon spectrum of both mono- 

(D . 

and bi-layer graphene systems. The results presented are also compared with those obtained in 

conventional 2DEG. The chiral nature of the quasiparticles in mono and bilayer graphene systems 

results in the observation of vr and 27r Berry's phase in the SdH- type oscillations in the plasmon 
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I. I. INTRODUCTION 

Recent progress in the experimental realization of both monolayer and bilayer graphene 
has led to extensive exploration of the electronic properties in these systems|l|, |2|. Ex- 
perimental and theoretical studies have shown that the nature of quasiparticles in these 
two-dimensional systems is very different from those of the conventional two-dimensional 
electron gas (2DEG) systems realized in semiconductor heterostructures. Graphene has a 
honeycomb lattice of carbon atoms. The quasiparticles in monolayer graphene have a band 
structure in which electron and hole bands touch at two points in the Brillouin zone. At 
these Dirac points the quasiparticles obey the massless Dirac equation leading to a linear 
dispersion relation e^ = vpk (with the Fermi speed vp = lO^m/s). This difference in the 
nature of the quasiparticles in monolayer graphene from conventional 2DEG has given rise 
to a host of new and unusual phenomena such as the anamolous quantum Hall effects and a 



m 



TT Berry phase[l|, |2|. These transport experiments have shown results in agreement with the 
presence of Dirac fermions. The 2D Dirac-like spectrum was confirmed recently by cyclotron 
resonance measurements and also by angle resolved photoelectron spectroscopy (ARPES) 
measurements in monolayer graphene [3|. Recent theoretical work on graphene multilayers 
has also shown the existence of Dirac electrons with a linear energy spectrum in mono- 
layer graphenejj]. On the other hand, experimental and theoretical results have shown that 
quasiparticles in bilayer graphene exhibit a parabolic dispersion relation and they can not 
be treated as massless but have a finite mass. In addition, The quasiparticles in both the 
graphene systems are chiral 2, m, la, la, lZ| . 

Plasmons are a very general phenomena and have been studied extensively in a wide 
variety of systems including ionized gases, simple metals and semiconductor 2DEG systems. 
In a 2DEG, these collective excitations are induced by the electron-electron interactions. 
Collective excitations (plasmons) are among the most important electronic properties of a 
system. In the presence of an external magnetic field, these collective excitations are known 
as magnetoplasmons. Magnetic oscillations of the plasmon frequency occur in a magnetic 
field. Single particle magneto-oscillatory phenomena such as the Shubnikov-de Haas and 
de Haas- van Alphen effects have provided very important probes of the electronic structure 



of solids. Their col 
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ective analog yields important insights into collective phenomena 
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Collective excitations of Dirac electrons in monolayer and bilayer 



graphene in the absence of a magnetic field have been investigated [l6|, llTI, [l8|, [l9|, |20| . 
Magnetic field effects on the plasmon spectrum have not been studied so far. In addition, 
since the quasiparticles in graphene are chiral, the particles will acquire Berry's phase as they 
move in the magnetic field leading to observable effects on the plasmon spectrum. To this 
end, in the present work, we study the magnetoplasmon spectrum within the self-consistent- 
field approach for both the monolayer and bilayer graphene systems. Magnetoplasmons can 
be observed by in-elastic light scattering experiments as revealed in studies carried out 
on 2DEG systems ll|, [l^, [l3|, [ij, [l5|. Similarly, in-elastic light scattering experiments 
are expected to yield information about the magnetplasmons in graphene. Furthermore, 
the results presented here can also be experimentally observed by Electron Energy Loss 
Spectroscopy (EELS) on graphene [21 1. 



II. ELECTRON ENERGY SPECTRUM IN MONOLAYER GRAPHENE 



We consider Dirac electrons in graphene moving in the x — y-plane. The magnetic field 
(B) is applied along the z-direction perpendicular to the graphene plane. We employ the 
Landau gauge and write the vector potential as A = (0, Bx, 0). The two-dimensional Dirac 
like Hamiltonian for single electron in the Landau gauge is (^ = c = 1 here) [l|, |2| 



Ho = VFO.{-iV + eA). 



(1) 



Here a = {cr^^, cTyjare the Pauli matrices and vp characterizes the electron Fermi velocity. 
The energy eigenf unctions are given by 



^n,fc,(r) 
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I = \J\/eB is the magnetic length, xq = I ky, Ly is the y-dimension of the graphene layer 
and Hn{x) are the Hermite polynomials. The energy eigenvalues are 

(3) 



e{n) = Ugy/n 



where Ug = v\f2eB is the cyclotron frequency of the monolayer graphene and n is an integer. 
Note that the Landau level spectrum for Dirac electrons is significantly different from the 



spectrum for electrons in conventional 2DEG which is given as e{n) = hujc{n + 1/2). The 
Landau level spectrum in graphene has ^/n dependence on the Landau level index as against 
linear dependence in 2DEG. The monolayer graphene has four fold degenerate (spin and 
valley) states with the n = level having energy e{n = 0) = 0. The quasiparticles in this 
system are chiral exhibiting vr Berry's phase. 



III. ELECTRON ENERGY SPECTRUM FOR BILAYER GRAPHENE 

The Landau level energy eigenvalues and eigenfunctions are given by 

e{n) = ujh\/n{ 
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where ± assigned to electron and hole states, tUfe = -^ is the cyclotron frequency of electrons 
in bilayer graphene and m* is the effective mass given as 0.044?7ie with rrie being the bare 
electron mass. The Landau level spectrum of electrons given by Eq.(4) is distinctly different 
from that of monolayer graphene and conventional 2DEG system. The electrons in bilayer 
are quasiparticles that exhibit parabolic dispersion with a smaller effective mass than the 
standard electrons. Bilayer graphene has four fold degenerate (spin and valley) states other 
than the n = level with energy e{n = 0) = which is eight-fold degenerate. These 
quasiparticles are chiral exhibiting 27r Berry's phase. 



A. INTER-LANDAU-BAND PLASMON SPECTRUM OF MONOLAYER AND 
BILAYER GRAPHENE IN A MAGNETIC FIELD 



The dynamic and static response properties of an electron system are all embodied in the 
structure of the density- density correlation function. We employ the Ehrenreich- Cohen self- 



consistent-field (SCF) approach [22(] to calculate the density-density correlation function. 
The SCF treatment presented here is by its nature a high density approximation which has 
been successfully employed in the study of collective excitations in low-dimensional systems 
both with and without an applied magnetic field. It has been found that SCF predictions 
of plasmon spectra are in excellent agreement with experimental results. Following the SCF 
approach, one can express the dielectric function as 

e(g, u) = l- Vciq)U{q, lu). (7) 

where the two-dimensional Fourier transform of the Coulomb potential Vc(q) = ^^, q = 
{ql+qyY^'^, K is the background dielectric constant and n(g, cu) is the non-interacting density- 
density correlation function 

X [5(n) - e{n) +uj + ir]]~^, (8) 

where C„„' (x) = (n2!/ni!) (x)"^ "^ e~^ L^^ "^(^) with rii = max(n, n'),n2 = min(n, n'), 

I -2 

L^{x) an associated Laguerre polynomial with x = ^^ here. This is a convenient form of 
n(g, uj) that facilitates writing of the real and imaginary parts of the correlation function. 
The plasmon modes are determined from the roots of the longitudinal dispersion relation 

l-Vc{q)ReU{q,uj) = (9) 

along with the condition Iinll{q,uj) = to ensure long-lived excitations. Employing Eq.(8), 
Eq.(9) can be expressed as 

and factor of 2 due to valley degeneracy. The plasmon modes originate from two kinds 
of electronic transitions: those involving different Landau bands (inter-Landau band plas- 
mons) and those within a single Landau-band (intra-Landau band plasmons). Inter-Landau 
band plasmons involve the local 2D magnetoplasma mode and the Bernstein-like plasma 
resonances, all of which involve excitation frequencies greater than the Landau-band sepa- 
ration. Since, in this work, we are not considering Landau level broadening therefore only 
the inter-Landau band plasmons will be investigated. 



We now examine the inter-Landau-band transitions. In this case n 7^ n'and Eq.(ll) yields 

/(g(n)) 

where A = {e{n) — s{n')) which permits us to write the following term in Eq (10) as 

(T( \^T( W o ^/(^H) (-^n^ 

{h{u) + hi-u)) = 2 ^^p_^^p - (13) 

Next, we consider the coefficient C„„'(x) in Eq.(lO) and expand it to lowest order in its 
argument (low wave- number expansion). In this case, we are only considering the n' = n±l 
terms. The inter-Landau band plasmon modes under consideration arise from neighboring 
Landau bands. Hence for n' = n + 1 and x <^ 1, using the following associated Laguerre 

n 

polynomial expansion L^{x) = y ^ (— 1)™ (i+rn)\(n-~m^\ 7^ for / > [23i] and retaining the first 

m=0 

term in the expansion for x <^ 1, C„„/(x) reduces to 

Cn,n+i{x) ^ {n + l)x, (14) 

and for n' = n — 1 and a; ^ 1, it reduces to 

Cn,n-i{x) -^ nx. (15) 
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Substitution of equations (13) and (14, 15) into equation (10) and replacing x = -^ yields 
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(16) 
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In obtaining the above result we note that A = {yrH — ^Jr\}j ujg. Therefore, A = -^j^ for 
n' = n + 1, and A = —y/e for n' = n — 1. We are considering the weak magnetic field 



case where many Landau levels are filled. In that case, we may substitute y/np for ^/n in 
Equation (16). '^f = ( ff ) is the Landau level index corresponding to the Fermi energy 
Ep. Equation (16) can be expressed as 



In terms of the 2D electron density n2D = ^^^/(^n) the inter-Landau-band plasmon 

n 

dispersion relation for monolayer graphene can be expressed as 

Corresponding calculation for bilayer graphene can be carried out. The equation that 
replaces Eq.(16), given above for monolayer graphene, is 

^ 27re2 2 ^/, , ,,, q\ 2{uJt)f{en) 



n ^ 
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For bilayer graphene Eq.(19) can be expressed as 



4:716 1 / m*ujhsr^ ,., ^ \ /^^v 



If we define n2D = "^^'' / J{^n) and the plasma frequency 



as 



2 A'Kn2De'^ _ , . 

^P,2D = ^9, (21) 

^' Kin* 

then the inter-Landau-band plasmon dispersion relation for bilayer graphene is 

u^={uj,f + ujl^o. (22) 



B. DISCUSSION OF RESULTS 

Eqs.(18) and (22) are the central results of this work. Eq.(18) is the inter-Landau band 
plasmon dispersion relation for monolayer graphene. The inter-Landau band plasmon energy 
as a function of the inverse magnetic field for the monolayer and bilayer graphene system 
with the plasmon energy for 2DEG at zero temperature is presented in Figs. (1,2). The 
following parameters were employed for doped graphene (Si02 substrate): k = 2.5, n2D = 
3 X 10^^ m^^ vf = lO^m/s. For the conventional 2DEG (a 2DEG at the GaAs-AlGaAs 
heteroj unction) we use the following parameters: m = .07me{Tne is the electron mass), 
K = 12 and n2D = 3 x 10^^ m"^. For electron density and magnetic field considered, electrons 
fill approximately 30 Landau levels, the upper limit in the summation for n2D is taken to be 
n = 30 while the lower limit is n = 0. In Fig.(l) we have plotted the plasmon energy as a 



function of the inverse magnetic field for both monolayer graphene and conventional 2DEG. 
The SdH-type oscillations are clearly visible that are a result of emptying out of electrons 
from successive Landau levels when they pass through the Fermi level as the magnetic field 
is increased. The amplitude of these oscillations is a monotonic function of the magnetic 
field. These oscillations have a tt Berry's phase due to the chiral nature of the quasiparticles 
in this system, the phase acquired by Dirac electrons in the presence of a magnetic field [l|. 
We also observe that the plasmon energy is ~ 4.2 times greater than in the 2DEG for the 
parameters considered. This is essentially due to the higher Fermi energy of the electrons 
in graphene and the smaller background dielectric constant. 

For bilayer graphene, we consider Eq.(22). There are two main differences between the 
plasmon dispersion relation for bilayer graphene given in Eq.(22) and the standard 2DEG 
result. Firstly, the cyclotron frequency Ub in bilayer is ~ 2 greater than the cyclotron 
frequency Uc at the same magnetic field in 2DEG due to the difference in the effective 
masses of the electrons in the two systems. Secondly, the 2D plasma frequency ujp^2D is 
also larger than in 2DEG for the same wave number q due to the smaller effective mass of 
electrons in bilayer compared to 2DEG and the smaller background dielectric constant k = 3 
in bilayer. The inter-Landau band plasmon energy as a function of the inverse magnetic 
field for doped bilayer and the 2DEG is shown in Fig. (2). The following parameters were 
used (Si02 substrate): k = 3, n2D = 3 x 10^^ m~^ and m* = 0.044meWith rrie being 
the usual electron mass. We again observe the SdH-type oscillations whose amplitude is a 
monotonic function of the magnetic field. We observe that the plasmon energy is ~ 2.6 times 
greater than in the 2DEG due to the smaller effective mass, valley degeneracy and smaller 
background dielectric constant. Due to the chiral nature of the quasiparticles in bilayer 
graphene, 2n Berry's phase is evident in the SdH type oscillations displayed in Fig. (2). 

In conclusion, we have determined the inter-Landau band plasmon frequency for both 
monolayer and bilayer graphene employing the SCF approach. The inter-Landau band 
plasmon energy is presented as a function of the inverse magnetic field. The SdH-type oscil- 
lations are clearly visible in both the systems and their amplitude is a monotonic function 
of the magnetic field. Due to the chiral nature of the quasiparticles in the mono and bilayer 
graphene system, vr and 2tt Berry's phases are observed in the SdH- type oscillations in the 
plasmon spectrum. 
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